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Soft processes in pp and AA interactions are considered in the framework of phenomenological 
model with color strings formation and fusion. Elementary parton collisions are realized in the 
model as the interaction of two colour dipoles from projectile and target nucleons. Modeling of 
the exclusive distributions of parton momentum fractions and transverse coordinates is performed. 
The interaction of colour strings in transverse plane is carried out in the framework of local string 
fusion model with the introduction of the lattice in the impact parameter plane and taking into 
account the finite rapidity length of strings. The parameters were fixed with experimental data on 
pp total inelastic cross section and charged multiplicity. 

The model was used for the calculation of long-range correlations between the multiplicities (n) 
and the mean transverse momenta (pt) of charged particles. The dependence of n-n, pt-n, pt-pt 
correlations on the width and position of the backward and forward rapidity windows was studied. 
Note that the model enables to describe the AA interactions without referring to the Glauber 
picture based on the concept of elementary nucleon-nucleon collisions. In this connection the 
charged multiplicity, the mean numbers of participant nucleons and binary collisions and their 
variances in the case of PbPb collisions were calculated and compared with the predictions of 
alternative models and the experimental data. The influence of different ways of centrality deter- 
mination on the multiplicity fluctuations and long-range correlations was also discussed. 
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1. Introduction 

The present work is devoted to development of Monte Carlo model for the description of soft 
proton-proton and heavy ion collisions at high energy. Due to inapplicability of usual QCD meth- 
ods in this region, such studies are semi-phenomenological. One of the models, which successfully 
describe particle production in soft region, is so-called quark-gluon string model [1,2]. The model 
involves two-stage scenario of particles production: at the first stage extended in rapidity objects 
(quark-gluon strings, or colour tubes) are stretched between the partons of colliding hadrons, and 
they fragment into observable particles at the second stage. 

In order to study these extended in rapidity strings it was proposed to measure long-range 
conelations between observables from two rapidity windows, separated by a gap. 

As these objects are finite in transverse plane, the overlapping of them and possible interactions 
can lead to non-trivial effects. Such interaction of the strings is considered in so-called string fusion 
model [3-5]. In particular, this model predicts the existence of non-zero long range correlations 
between multiplicity and mean transverse momentum (pt-n) and pt-pt correlations. The theoretical 
study of long-range correlations in the framework of string fusion model is carried out using Monte 
Carlo simulations [6- 10]. The goal of this work is to develop detailed Monte Carlo model, taking 
into account finite rapidity width of the strings and their fusion at the transverse plane both for pp 
and AA collisions. 

2. The model 

In order to provide parton-string model, which is suitable for reliable description of fluctua- 
tions and correlations, one should take into account not only inclusive, but also exclusive parton 
distributions as well as the details of the string formation and fusion. 

2.1 Partonic picture of p-p collision 

The form of inclusive distributions on the momentum fraction for N=2n partons is taken from 



At « > 1 the sea quarks and antiquarks have the same distribution as the valence quarks. Poisson 
distribution for the number of quark-antiquark (diquark) pairs {n) is assumed with some parameter 
A, omiting the case of « = 0, 

For corresponding exclusive distribution for an arbitrary number of quark-antiquark pairs we 



[11, 12]: 
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find: 
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The valence quark is labelled by N- 1 , the diquark - by N, and the other refers to sea quarks 
and antiquarks. = 3/2 (ud-diquark), aiM = 5/2 (uu-diquark). 

Corresponding algorithms for modelling of such distributions are discussed in [10]. 

The exclusive distribution in the impact parameter plane is constructed on the following sup- 
positions: 

N 

1 . The position of the mass centre is fixed: X! " = 0- 

2. The inclusive distribution of each parton is the 2-dimensional Gaussian distribution. 

N 

3. The normalization condition is < >=< L ''/^ >= ''o^, 

where the parameter ro^ is connected with the mean square radius of the proton: < >= ^Tq. 

For the description of inelastic proton-proton collisions we assume that the elementary parton 
collision is implemented as an interaction of two color dipoles consisting of a valence quark and 
diquark, or of a quark-antiquark pair. The probability amplitude of interaction of two dipoles with 
coordinates (rir2) and (r^r^) in the impact parameter plane is given by [13]: 

where as - is a. constant. 

The value of as is assumed to be an effective coupling constant, and as is used as a fitting 
parameter for better description of the experimental data, it is assumed that the value does not 
depend neither on the energy, nor on the number of quark-antiquark pairs in the proton. 

Taking into account the confinement effects [13, 14] leads: 



■' 2 



^0 I ]+kJ ]-kJ -^n ^ "'^ - '"^1 



(2.3) 



which provides an exponential decrease of collision probability of two protons with large impact 
parameter. 

In the eikonal approximation [14, 15] the probability of interaction of two dipoles is given by: 



Pij 



l-e^f". (2.4) 



-E/.7 

The total probability of inelastic interaction of two ptorons is p = l—e , where the summation 
is made over all the dipoles. 

For more detiled discussion please refer also to [9, 10]. 

2.2 Calculation of observables 

Probabilities of dipole interaction, obtained in the previous section, are used for the construc- 
tion of the collision matrix. It should be noted, that a dipole can interact only with one other dipole, 
so filling of the interaction matrix starts from the valence dipole. 

The next step is the generation of strings in rapidity space. Rapidity ends of a string ymm,ymax 
are determined from the kinematic condition of a string decay only on two particles with an average 
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Pt = 03GeV and masses = O.lSGeV (for pion - in case of quarl-antiquatk string) or nip = 
0.94GeV (for nucleon - in case of diquark at the end of the string). 

The strings, that are too short, are excluded from consideration by the requirement that the 
sum of the masses of the particles produced should be less than the mass of the string, ie ^/sxAXB, 
where xa,xb - momentum fractions of the partons at the ends of the string. Transverse coordinates 
of the center of the strings are equal to the arithmetic mean of the corresponding coordinates of the 
partons at the ends. 

Due to finite transverse size of the strings they overlap in the impact parameter plane. The 
interaction of colour strings in transverse plane is carried out in the framework of local string 
fusion model [5] with the introduction of the lattice in the impact parameter plane. [6, 7, 16- 18]. 
Cellular variant of string fusion has been chosen in order to simplify taking into account finite 
rapidity length of strings, because, as it had been shown, numerical results differ little in several 
string fusion model variants [6, 8]. 

According to this model, mean multiplicity of charged particles and mean pt originated from 
the area S^, where k strings are overlapping are the following: 

(Ai), = Mo\^— {P,)k = Po<^, (2.5) 
Co 

where Gq = Tlr^^ ~ string transverse area 

In the discrete model transverse plane is considered as a grid with the cell area equals to string 
transverse area and strings are fused if their centers occupy the same cell. 

Finite rapidity width of strings is also taken into account[10]: rapidity bins are defined by 
strings edges, so every rapidity interval contains constant number of strings, then bins are processed 
separately with the final summation over whole rapidity axis. 

The Monte Carlo model was implemented as C++ class, and in order to get results with high 
statistics, distributed grid computing have been involved [19]. 

Parameters fixation 

In the framework of our model the growth of multiplicity and total inelastic cross section is 
achieved by increasing the number of quark-antiquark pairs (parameter A). The value of A is turned 
from the experimental data on the total inelastic cross section [20]. 

The simulations are conducted in two phases. The first step is devoted to determination of 
(7pp for the given A. Using experimental data, the dependence of X on collision energy is obtained 
(example is shown at fig. 1). 

One should note the importance of sea quark presence even at low energy [9] . 
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10' 'O' ■is, GeV 



Figure 1: Connection of parameter A and the energy for the case of ttj = 0.7. 

At the second step all other calculations are carried out for the given energy. 
The value of mean multiplicity per rapidity from one single string jUq is fixed at one point at 
low energy. 

The other parameters of the model were chosen as follows: for rmax — 0.2 — 0.3 fm; ra- 
tio rmax/ro is chosen to be 0.5; constant as is rixed for best description of the dependence of 
charged multiplicity on energy [21]. For the string radius r^^r we considered several possible val- 
ues: 0.2/m,0.3/m,0.4/m, which are within the estimates [22, 23, 16, 17], used for the description 
of long-range coiTclations in nucleus-nucleus collisions. 

The results of parameters fixing are presented in Table 1 , the corresponding plot for the multi- 
plicity shown in Fig. 2. 



Table 1: 



rstr=0.2fm 


a5=0.7 


Ato=0.92 


rstr=03fm 




Ato=1.02 


r^,^=0.4fm 




Ato=l-12 
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A ISR (pp) INEL 
• UA5{pp)INEL 
T PHOBOS (pp) INEL 
« CDF (pp) NSD 
■ ALICE (pp) INEL 

— MC Model (rstr=0-2 fm. QS=0.7) 

— MC Model (rstr=0.3fin. as=0.5) 
— - MC Model (rstr=0.4 fm, QS=0.5) 



UA1 (pp) NSD 
UA5 (pp) NSD 
STAR (pp) NSD 
CMS (pp) NSD 
ALICE (pp) NSD 




I 1 1 1 
10^ 



, ^ 

Energy n/^ (GeV) 



Figure 2: Chaged-particle multiplicity in p-p collisions. Calculations of the Monte Carlo model and exper- 
imental data [21,24]. 



The calculation results show good description of the multiplicity in a wide energy range, and 
taking into account the results of the LHC at 7 TeV we use the first set of parameters {rgtr = 0.2fm, 
as = 0.7, jiQ = 0.92) for further studies. 

Calculation of correlation functions 

The correlation function between two observable B (backward) and F (forward) from different 
rapidity windows is defined as the mean value of B with fixed F [6] : 

fsMF) = {B)F. (2.6) 
The correlation coefficient represents the slope of correlation function: 

bs-F = \f=<f>- (2.7) 

dF 

Often it is useful to switch to normalized variables: B^B/ {B),F F/{F), in this case pt — n 
correlation coefficient would become dimensionless, and both n — n and pt — pt correlation coeffi- 
cients would not change in case of symmetrical windows. Thus we use the following definitions: 

_ <nF> d<nB> 
< ng > drif 



<nf > d < ptB> : 



(2.9) 



^p,-p, - / „ ^ ■ \prf=<p,F>- (2-10) 



Pt-n ^ „ ^ j„ \nF=<nf>i 

< ptB > drif 

< PtF > d < PtB > 

< PtB > dptp 

In practice, the experimental determination of the correlation coefficient and Monte Carlo 
simulation is performed by obtaining fitting of f{F) by linear function from (F) — a^^ ^ to {F) + 

where a^^ ^ = -J < F^ > — {F)^. This method is used in this paper. 
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For the calculation of these correlations in the framework of string fusion model the occupancy 

T], of each cell is defined from string configuration and one assumes, that every cell of the lattice 

emit particles according to Poisson distribution with parameter equal to T], • pto^y. The mean event 

transverse momentum is generated with Gaussian distribution with mean [7] (£«iO;i)/(£«,) and 

(' (' 

variance Y^niDi, where a, = i^/po, D, = y^D, D = jpo. 

i 

The coefficient 7 is the coefficient of proportionality between the average transverse momen- 
tum and root of the variance of transverse momentum from one single string. This parameter is 
assumed to be independent of the energy. As in [7], we used the value 7 = 0.61. 

Obtained by this way mean event transverse momentum and multiplicity in both forward and 
backward windows ai^e used to will 2d histogram (fig. 14) and correlation function is obtained by 
regression. 

This way in case of n-n and pt-n correlations is fully equivalent to the one, described in [6, 10] 
and also coincides with the one, used for pt-pt congelations [7]. 

3. Results on pp interaction 
Correlation functions 

Correlation functions for p-p collisions are sown at fig. 3, 4 and 5. One should note the non- 
Unearity of all types correlation functions with significant deviation from linear function occuring 
only outside the region {np) — Of <np < (np) + af, which is used for determination of correlation 
coefficients. 



n-n correlation function 




<n^> 



Figure 3: n-n correlation function for pp collisions at 7 TeV, calculated in MC model. Rapidity windows 
are (-0.8, 0) and (0 ,0.8). 
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p^-n correlation function (-0.8, 0), (0, 0.8) 




Figure 4: pt-n correlation functions for pp collisions at 7 TeV, calculated in MC model. 
Pj-Pj correlation function (-0.8, 0) (0, 0.8) 




<ptp> 

Figure 5: pt-pt correlation functions for pp collisions at 7 TeV, calculated in MC model. 

Correlation coefficients 

We also studied the dependence of correlation coefficients on the width of rapidity windows. 

Firstly we consider the dependence of correlation coefficients on the width of backward rapid- 
ity window with fixed forward window at (0.6, 0.8). The configurations of the windows used are 
shown at the fig. 6. Correlation coefficients are calculated for pp collisions at 7 TeV and results are 
shown at the fig. 7. 
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0.6 0.8 



Figure 6: Configuration of rapidity windows with fixed forward window. 
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Figure 7: n-n, pt-n and pt-pt correlation coefficients in p-p collisions at 7 TeV as a function of the width of 
backward window. 

The dependence of correlation coefficients on the width forward rapidity windows is shown 
at fig. 9. 
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-0.8 



-0.6 



Figure 8: Configurations of rapidity windows with fixed backward window. 




Figure 9: n-n, pt-n and pt-pt correlation coefficients in p-p collisions at 7 TeV as a function of the width of 
forward window. 



The values calculated in the Monte Carlo program are fitted by the formula 
Fit parameters, obtained for n-n, pt-n and pt-pt correlations are shown in the table 2. 



(3.1) 
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Table 2: Parameters of a fitirovaniye of coefficients of correlation 



Type of correlation 




k 


n — n 


1.000 ±0.002 


0.34 ±0.01 


Pt-n 


0.088 ±0.002 


0.41 ±0.02 


Pt-Pt 


0.53 ±0.05 


4.6 ±0.6 



Discussion of results 

The absence of dependence of n — n correlation coefficient on the width of a backward win- 
dow corresponds to the predictions of the models of independent emitters [25] and to asimptotical 
expressions, obtained in the string fusion model [6, 7, 26] in the limit case of high string density. 

Small dependence of pt —n correlation coefficient on the width of a forward window can be 
caused by dependence of this coefficient on the gap between windows (which is stronger for pt — n 
correlation, than for n — n). An additional source of such dependence may be constraint «b > 0. In 
fact, if we select an event with non-zero number of particles in backward window, we potentially 
select higher-multiplicity event, probably, with higher mean transverse momentum and this can 
change the value of the correlation coefficient in the relative variables. In wide backward window 
this constraint becomes negligible. 

Pt — Pi correlations also do not depend on width of a back window within error bars. 

The dependence of a coefficient of correlation on the width of a forwai^d window is completely 
described by the formula (3.1); for n — n correlation j3 = 1 is found, that is in full agreement with 
the predictions of the model of the independent emitters [25] . The fact, that the value k for n — n 
and Pt — n coiTclations does not coincide, as it was predicted for high string density [6], means that 
the density of strings, where the asymptotic formula is applicable, is not reached in p-p collisions 
at 7 TeV. 

One should notice that pt — Pt correlations strongly depend on width of forwai^d rapidity win- 
dow, and at wide sufficiently Ayp dominate over pt — n correlations. 

From the physical point of view that fact that correlation coefficients depend only on the width 
of a forward window means [16, 23] that the dynamic variable in the forward window is the one 
that classifies events, and the wider window is, the closer this classification coincides with the one 
by string configurations, and hence, the more this value is correlated with a variable in another 
rapidity window. 

In addition, one should notice that obtained dependence of correlation coefficients on the width 
of rapidity windows can be confirmed in experiment; this is done for n — n correlations [27, 28], 
and used for the efficiency corrections. 
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Dependence of correlation coefficients on the gap between rapidity windows 

The dependence of correlation coefficients on the provision of rapidity windows was also 
studied in the model. The symmetric windows of the width 0.8 rapidity units were chosen. The 
results are shown at the fig. 10. 



-Q 0. 





ygap 



Figure 10: Dependence of n-n, pt-n and pt-pt correlation coefficients on the gap between rapidity windows 
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All three types of correlations show similar behavior and decrease at high gap in the same way, 
as rapidity distribution of charged particles. 

It is noticeable that both pt — n and pt — pt correlation coefficients reach zero at wide enough 
rapidity gap. It can be it is explained as follows: only the lengthiest valence strings contribute to 
these windows far away from midrapidity, therefore there is no strings to fuse in this area, and 
consequently there are no pt — n and pt — pt correlation. Multiplicity-multiplicity correlations 
remain due to fluctuation of the numbers of strings (volume fluctuation): zero and one. 

4. Nucleus-nucleus collisions 

For the description of nucleus-nucleus collisions at high energy the Glauber model is is widely 
used [29, 30], which is based on an assumption of incoherent superposition of single nucleon- 
nucleon collisions. It is supposed that trajectories of nucleons in a nuclei are approximately straight 
lines and all consequent collisions are occurred at the same cross section and with identical mean 
number of produced charged particles (the same as for pp collision). 

Thus, multiplicities in pp and AA collisions are connected by the following: 

N^=NcoiNPP. (4.1) 

In such manner Glauber model considerably overestimate charged multiplicity in comparison 
with experimental data and contains obvious breaking of the energy conservation: indeed, the in 
such picture same part of energy can go to the particle production several times [31]. 

In order to get an agreement with experimental data one should use instead of (4.1) some 
interpolations: 

A^,;, = C{xNp,rt + ( 1 - x)Ncon ) . (4.2) 

With turned parameters C and x one can achieve of distribution of charged multiplicity consistent 
with experimental data. The Glauber model is used in such manner for the definition of number of 
participating nucleons in those experiments where it is impossible to determine it directly [32]. 

This problem is considered also in modified Glauber model [31], where energy loss is taken 
into account effectively. There are also several models [33-35] representing "The Glauber on a 
parton level". 

In the present paper we develop a Monte Carlo model without refeiTing to the Glauber pic- 
ture based on the concept of elementary sequential nucleon-nucleon or partonic collisions. We 
implement direct generalization of pp model, described above. 

The initial aiTangement of nucleons is done with standard Woods-Saxon distribution: 
PA(r) = i+,,pg_^)/^| wifli/? = 1.07 -A'/'frn, d = 0.545/m. 

Nucleons are treated as set of the dipoles; the elementary collisions it is carried out by means 
of formulas (2.3 - 2.4). The nucleon is considered as participant if at least one of its dipoles faced 
a dipole of other nucleus. There is no additional parameter turning. The set of observables is 
calculated in the same way as for proton-proton collisions. 
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Figure 11: Number of participant nucleons, number of binary collisions (left), their variances (middle), 
scaled variances and correlator 2{{N pan part b) ~ {NpartA){NpartB)) (right): our Monte-Carlo model without 
Glauber picture (up), and calculations in the Glauber model at Of^i^ = 34mb (down), 
note: additional scale (at the right) is used for variables, marked by star. 

In the framework of the MC model we calculated the number of participants, binary collisions 
and their variances for PbPb collisions at low energy (Gnn = 34mb) at fixed impact parameter and 
compare them with results of Glauber model calculations (fig. 11). 

The results show practically the complete coincidence of predictions of two models on the 
number of participated nucleons, while the number of binary collisions is almost twice less. There 
is qualitatively similar behavior of variances of number of participants and number of collisions in 
these two models; quantitatively in the Glauber model fluctuations appear higher. 

The calculations carried out by the author in the Glauber model are close with analytical and 
numerical calculations [36, 37] Onn = 3l.4mb. 

The decrease of the number of binary collisions can be explained by the following: in every 
nucleus-nucleus collision each dipole can interact with other one only once, and the energy of 
colliding partons directly goes to particle production. It reduces additional probability of further 
interactions of the nucleons with akeady "used" partons. Thus, such picture respects the energy 
conservation. 

The charged multiplicity per rapidity over a half of number of participating nucleons at the 
LHC energy is shown at the fig. 1 2 and compared with ALICE experimental data. 

In general, we achieved a good agreement, remaining discrepancy probably could be reduced 
by additional parameters turning, taking into account some results on p-A or AA scattering. 
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•o Pb-Pb 2.76 TeV ALICE H MC model 
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Figure 12: Dependence of {dNc-i,/dri)/{Npan/'2) on the number of participants for Pb-Pb collisions at 
^s^VA' = 2.76TeV. Results of MC calculations and experimental data [32]. 



On fig. 13 long-range correlation coefficients of n — n, pt — n and pt — Pt correlations are shown 
for PbPb collisions at the energy of LHC (2.76 TeV). Two different configurations of forward and 
backward rapidity windows were considered: windows in width 0.8 and 2 rapidity units. 

It is noticeable, that even at the fixed impact parameter n — n correlation are very strong. The 
correlation coefficient decreases in the central collisions that it expected due to string fusion effects 
(decrease of n-n correlation) are important in central PbPb collisions. High value of n-n correlations 
in semi-central and peripheral collisions are due to volume fluctuation effects. 

The large value of pt — Pt correlations is remarkable and dominates over pt — n correlations. 
This fact is in agreement with predictions [7] stated that at LHC energies pt — pt correlation coef- 
ficient can reach the level of n — n correlations. 
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5. Conclusion 

Long-range correlations are calculated in the framework of the string fusion model, taking into 
account finite rapidity length of the strings. 

N-n, pt-n and pt-pt correlation functions were calculated and it was shown that they demon- 
strate non-linear behaviour. The dependence of the correlation coefficients on the width of rapidity 
windows and the gap between them is calculated and compared with the predictions of the model 
of independent emitters. 

The developed model enabled to describe also AA interactions without referring to the Glauber 
picture of nucleon-nucleon collisions. The number of soft binary collisions in this approach is 
proved to be less than in the Glauber approach. 

Correlation coefficients for Pb-Pb collisions are estimated in case of fixed impact parameter. 
In Pb-Pb collisions at LHC energy the strength of pt-pt correlation is larger compared to pt-n 
correlation. 
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Appendix 

At fig. 14 an example of n-n correlation cloud is shown. The colour denotes the number of 
events with given np and hb- N-n correlation function is obtained by regression of this histogram. 

n-n cloud 




Figure 14: 2d n-n correlation histogram, calculated in Monte Carlo model for pp interactions at 7TeV, 
rapidity windows are (-0.8, 0), (0, 0.8). 
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